AD-A126  123  WAVE  PROPAGATION  IN  PARTICULATE  MEDIA:  MULTIPLE 

SCATTERING  OF  OPTICAL  PULSES  IN  SCALE  MODEL  CLOUDS(U) 
OREGON  GRADUATE  CENTER  BEAVERTON  R  A  ELLIOTT  08  MAR  83 
UNCLASSIFIED  NOOO14-79-C-0897  F/G  20/14 


6123 


1  M-*, 


Unclassified 

iRiTt  ClWlHOTWW  or  ImH  Mg|  gMw  P».»  KwtfC 

REPORT  DOCUMENTATION  PACE  I  •eforeDcompletiTkgNfoKM- 

IPORT  NUMBER  I*  0OVT  ACCESSION  NO,  *ECipTEnT*A  CATALOG  NUMRER 


-fliaUi  3l 


14  TITLE  r«ntf  WMNMt 

WAVE  PROPAGATION  IN  PARTICULATE  MEDIA? 
MULTIPLE  SCATTERING  OF  OPTICAL  PULSES  IN  SCALE 
MODEL  CLOUDS 


ft.  TYPE  or  RCP©AT  ft  PERtOD  COVERED 

Technical  Report 
9/1/79-8/31/82 
'ft.  PK9tro*MiM0  006.  REPORT  numRER 


T.  AuTmOR<r> 

Richard  A.  Elliott 

N00014-79-C-0897 

•  RCRTORMIMO  organization  NAME  and  address 

Oregon  Graduate  Center 

19600  N.W.  Walker  Road 

Beaverton,  Oregon  97006 

10.  PROGRAM  Element  PROJECT.  TAIK 
AREA  ft  WORK  unit  WUMRERft 

NR394-045 

II.  CONTROLLING  OFFlCC  NAME  AND  ADDRESS 

Office  of  Naval  Research,  Physics  Program  Office 

•1.  RERORT  DATE 

March  8,  1983 

800  N.  Quincy  Street 

Arlington,  VA  22217 

is  number  or  rages 

60 

1ft  MONITORING  agency  name  ft  AODRE ftftflf  Pl//#reni  Itmm  CertirefUaft  Othf) 

1ft  ftECuRiT  Y  CLAftft  f#f  tAie  report; 

Unclassified 

IftA.  DECk.AftftirtCATlON  DOWNGRADING 
SCHEDULE 

1ft  DiftTRiRuTiON  STATEMENT  fef  thi»  Report; 

Approved  for  public  release;  distribution  unlimited. 


17.  DlSTRiiuTlOW  STATEMENT  (of  CM  eftefrocr  w%ff*  in  flee*  20,  II  4ttl*rwtt  frmm  *•*•*) 


DTIC 

ELECTE 


1ft  KEY  WORDS  (CriIIiih*  RR  nvpm  «<!•  II  mcmimt  ftp  blmtO  ^ 

Optical  scattering,  pulse  delay,  pulse  stretching,  beam  spreading,  cloud 
simulation 


so  abstract  •«  »•*•»•••*<•  n *•<•••«?  •*«•*«*  »»  w»»*  The  results  of  experiments  t 

measure  the  delay,  temporal  stretching,  attenuation  and  spatial  spreading  of 
optical  pulses  in  scale  model  clouds  are  reported.  The  model  clouds  consisted 
of  diiodomethane/water  or  paraffin  oil/water  emulsions  maintained  in  a  rotat¬ 
ing  scattering  cell  to  prevent  settling  of  the  droplets.  The  optical  pulses 
were  532  nm,  25  ps  duration  pulses  generated  by  a  frequency  doubled,  passively 
mode  locked  Nd:YAG  laser  and  were  detected  with  a  10  ps  resolution  streak 
camera.  The  measurements  of  the  delay  in  the  mean  arrival  time  of  the  pulses 
due  to  multiple  scattering  are  the  first  ever  measured  directl 


DO  M7J  t^TiOROf  f  ROv*»i»os*oi.tT«  Unclassified  N  _ 

&  O  0  o  p  g  Q  ^^ICUlllTT  CkAUl'lCATlOM  or  TniS  R  ABC  *"*•'•* 


WAVE  PROPAGATION  IN  PARTICULATE  MEDIA: 


MULTIPLE  SCATTERING  OF  OPTICAL  PULSES  IN  SCALE  MODEL  CLOUDS 


by 


Richard  A.  Elliott 


Technical  Report 
March  8,  1983 


Contractor:  The  Oregon  Graduate  Center 

Sponsor:  Office  of  Naval  Research 

800  N.  Quincy  Street 
Arlington,  Virginia  22217 


Aoooasloa  For 

mis  GRAM 

DTIC  TAB 

Uhannauncsd 

Justification. 


□ 


Distribution/ 


Availability  Codes 
*”  jAvail  and/or 


Dist 


A 


Special 


Contract  Number:  N00014-79-C-00897 


Effective  Date  of  Contract:  September  1,  1979-August  31,  1982 
Principal  Investigator:  Richard  A.  Elliott 


Reproduction  in  whole  or  in  part  is 
the  United  States 


permitted  for  any  purpose  of 
Government . 


Approved  for  Public  Release; 
Distribution  Unlimited 


This  manuscript  has  been  submitted  for  publication 
to  Applied  Optics* 


MULTIPLE  SCATTERING  OF  OPTICAL  PULSES  IN  SCALE  MODEL  CLOUDS* 

by 

Richard  A.  Elliott 
Oregon  Graduate  Center 
19600  N.W.  Walker  Road 
Beaverton,  Oregon  97006 

ABSTRACT 

The  results  of  experiments  to  measure  the  delay,  temporal 
stretching,  attenuation  and  spatial  spreading  of  optical  pulses  in 
scale  model  clouds  are  reported.  The  model  clouds  consisted  of 
diiodomethane/water  or  paraffin  oil/water  emulsions  maintained  in 
a  rotating  scattering  cell  to  prevent  settling  of  the  droplets. 

The  optical  pulses  were  532  nm,  25  ps  duration  pulses  generated  by 
a  frequency  doubled,  passively  mode  locked  Nd:YAG  laser  and  were 
detected  with  a  10  ps  resolution  streak  camera.  The  measurements 
of  the  delay  in  the  mean  arrival  time  of  the  pulses  due  to 
multiple  scattering  are  the  first  ever  measured  directly. 

*Part8  of  this  work  were  reported  at  the  Annual  Meetings  of  the 
Optical  Society  of  America  at  Kissimmee,  Florida  October  29,  1981 
and  Tucson,  Arizona,  October  21,  1982;  at  the  Conference  on  Lasers 
and  Electro-Optics,  Phoenix,  Arizona,  April  16,  1982;  and  at  the 
OSA  Topical  Meeting  on  Optical  Techniques  for  Remote  Probing  of 
the  Atmosphere,  Incline  Village,  Nevada,  January  11,  1983. 


MULTIPLE  SCATTERING  OF  OPTICAL  PULSES  IN  SCALE  MODEL  CLOUDS 
I.  Introduction 

Multiple  scattering  by  aerosols,  dust,  and  water  droplets  in 
clouds  or  fogs  is  responsible  for  a  number  of  deleterious  effects 
which  severely  impact  the  performance  of  imaging,  remote  sensing, 
LIDAR  and  optical  communication  systems.  The  work  presented  here 
bears  some  relevance  for  other  systems,  but  is  mainly  concerned 
with  the  effects  of  multiple  scattering  on  laser  pulse 
communication  systems  such  as  the  blue-green  satellite  to 
submarine  system  of  interest  to  the  U.S.  Navy. 

The  past  decade  has  seen  a  considerable  effort expended 
on  understanding  optical  pulse  propagation  in  clouds  and  fog 
beginning  with  Bucher's  pioneering  experiments*  and  Monte  Carlo 
simulations^  and  extending  to  the  recent  experiments  on  the 
island  of  Kauai. This  work  has  shown  the  dramatic  effects 
which  are  possible  when  multiple  scattering  and  multiple  path 
effects  can  play  a  dominant  role  as  occurs  in  optically  thick 
clouds.  Nanosecond  duration  pulses  have  been  observed  to  be 
temporally  stretched  to  more  than  microsecond  duration  and 
theoretical  models  predict  a  comparable  time  delay.  An  initially 
collimated  laser  beam  may  be  spread  in  spatial  extent  by  an  amount 
approximat ing  the  geometrical  thickness  of  the  cloud  and  can  have 
a  nearly  uniform  angular  distribution.  The  only  effect  which 
could  be  said  not  to  be  deleterious  is  the  deviation  of  the 


attenuation  of  the  transmitted  pulse  from  the  Beer-Laabert  lav: 
use  of  a  vide  field  of  viev  receiver  may  allov  collection  of  as 
much  as  100  dB  more  energy  than  predicted  by  the  single  scatter 
theory.^ 

The  current  state  of  understanding  as  represented  by  existing 
theoretical  models  and  scaling  lavs  provides  a  basis  for  a  good 
description  of  the  effects  of  multiple  scattering.  The  available 
experimental  evidence  although  incomplete  does  agree  reasonably 
well  with  these  models  and  serves  to  lend  validity  to  them. 

A  common  characteristic  of  most  scattering  experiments 
performed  in  the  field  on  naturally  occurring  phenomena  is  the 
inability  to  control  experimental  conditions  and  the  difficulty 
involved  in  characterizing  the  scattering  medium.  It  is  not  in 
general  possible  to  adequately  instrument  the  entire  propagation 
path  to  monitor  variations  in  the  medium  nor  to  ensure  the 
reasonable  range  of  conditions  necessary  to  validate  a  theory. 
Laboratory  experiments  thus  provide  a  much  needed  supplement  to 
the  data  acquired  in  field  work. 

The  laboratory  system  which  should  logically  be  used  to  study 
scattering  in  clouds  is  an  artificially  generated  water  droplet 
cloud  since  this  provides  the  greatest  degree  of  similarity 
between  the  natural  and  model  systems.  Such  experiments 
have^-^  indeed  provided  very  useful  information  on  scatterir_ 
in  optically  thin  clouds  but  since  the  number  density  of  water 
droplets  which  can  be  maintained  even  under  the  best  conditions  is 


not  much  greater  than  that  in  atmospheric  clouds  and  the  overall 
dimensions  of  artificial  clouds  is  on  the  order  of  a  few  meters  an 


optical  thickness  of  10  is  near  the  maximum  obtainable.  In  order 
to  study  scattering  in  optically  thick  media  in  the  laboratory  it 
is  a  practical  necessity  to  employ  other  model  systems. 

A  model  system  should  be  stable  for  extended  periods  and  the 
important  scattering  parameters  subject  to  some  measure  of  control 
and  be  readily  characterized.  A  commonly  used  system  for 
laboratory  studies  of  Mie  scattering,  i.e.,  scattering  from 
particles  with  diameters  comparable  to  or  larger  than  the 
wavelength ,  is  a  suspension  of  latex  spheres  in  water. 20  These 
are  generally  available  in  monodisperse  size  distributions  with 
diameters  ranging  from  tenths  to  hundreds  of  micrometers.  The 
quantities  needed  for  multiple  scattering  experiments  involving 
optical  thicknesses  in  excess  of  100  are  however  prohibitively 
expensive  and  moreover  only  a  single  ;ndex  of  refraction  is 
accessible. 

An  alternative  to  latex  spheres  is  provided  by  emulsions  of 
dielectric  liquids.  For,  one  liquid  forms  spherical  droplets 
dispersed  throughout  the  other  and  the  optical  properties  of  these 
systems  depend  entirely  on  the  liquids  used  and  the  size 
distribution  of  the  dispersed  droplets.  The  size  distribution  is 
remarkably  stable  especially  when  the  emulsion  is  formed  in 
conjunction  with  a  surfactant  and  a  wide  range  of  refractive 
indices  is  available.  For  this  reason  all  the  experiments 


reported  here  were  performed  on  two  types  of  emulsions,  either 
paraffin  oil  or  diiodomethane  dispersed  in  water.  The  latter 
system  is  of  special  interest  since  the  relative  index  of 
refraction  of  diiodomethane  to  water  is  nearly  the  same  as  that  of 
water  to  air  and  hence  is  a  good  model  of  a  cloud. 

The  mathematical  models  describing  the  effects  of  multiple 
scattering  on  pulse  propagation  give  expressions  for  the  pulse 
stretching,  delay,  attenuation  and  beam  broadening  in  terms  of  the 
optical  thickness,  the  scattering  length  and  single  scatter 
parameters  of  the  medium.  It  is  possible  through  the  use  of  an 
appropriate  emulsion  to  duplicate  the  single  scatter  properties  of 
a  cloud.  The  optical  thickness  is  the  geometrical  thickness 
divided  by  the  scattering  length,  idiich  is  just  the  reciprocal  of 
the  number  density  times  the  average  scatter  cross  section.  Thus 
the  optical  thickness  is  directly  proportional  to  number  density 
and  use  of  a  dense  medium  makes  it  possible  to  reduce  the  physical 
dimensions  of  even  optically  thick  systems  to  laboratory  scale. 

The  scattering  length  or  optical  mean  free  path  is  however 
inversely  proportional  to  number  density.  In  a  dense  system  then 
the  time  interval  between  scattering  events  is  reduced  and  faster 
detectors  and  shorter  duration  incident  pulses  must  be  employed. 
For  example,  15  ns  duration  optical  pulses  are  stretched  to  ~  10 
Us  duration  in  clouds  a  few  km  thick. ^*^,6  To  replicate  this 
situation  on  a  laboratory  scale  the  scattering  length  must  be  at 


least  10  times  smaller  and  the  detection  system  must  be  10  times 
faster.  This  was  achieved  in  the  experiments  reported  here  by 
using  a  passively  mode-locked  Nd:YAG  laser  frequency  doubled  to 
532  nm  as  the  pulse  source  and  a  streak  camera  capable  of  10  ps 
time  resolution  as  the  detector. 

II.  Theoretical  Models 

Mathematical  models  of  the  effects  of  multiple  scattering  on 
pulse  propagation  have  been  derived  from  computer  simulations  and 
analytic  solutions  to  simplified  forms  of  the  equation  of 
radiative  transfer. 2,8-15  The  physical  situation  most 
frequently  modelled  is  that  of  a  cloud  layer  of  thickness  Z  in  the 
propagation  direction  and  infinite  in  extent  in  the  transverse 
plane . 

The  Monte  Carlo  simulations^  proceed  by  assuming  a  ray 
enters  the  top  of  the  layer,  proceeds  a  distance  encounters  a 
scattering  center  which  it  leaves  in  a  direction  specified  by 
angles  0i  and  4>i  relative  to  its  incident  direction,  proceeds  a 
distance  and  is  again  scattered.  At  each  scattering  a  new 
direction  (0j,  relative  to  its  direction  before  that 
scattering  is  specified.  The  distances  between  scatterings, 
are  distributed  exponentially;  the  azimuthal  angles  are 
uniformly  distributed  between  0  and  2n;  and  the  polar  angles 
are  distributed  according  to  the  single  scatter  function,  P(6), 
appropriate  for  the  size  distribution  and  index  of  refraction  of 
the  particles  of  the  system  being  modelled. 


The  ray  is  propagated  in  this  manner  until  it  exits  from  the 
'cloud'  boundary.  The  path  length,  £&£;  the  displacement  from 
the  axis  and  angle  of  exit  are  recorded  and  the  procedure  repeated 
through  many  trials.  It  is  then  possible  to  determine  the  mean 
and  variance  of  the  path  length,  and  the  spatial  and  angular 
distribution  of  the  exiting  rays.  Bucher^  found  that  the 
multipath  effects  in  thick  clouds  could  be  accurately  scaled 
according  to  simple  equations  involving  only  a  few  parameters:  the 
optical  thickness  of  the  cloud,  T;  the  mean  distance  between 
scatterings,  b;  and  the  asymmetry  factor,  g.  These  parameters  are 
defined  in  terms  of  the  geometrical  dimensions,  number  density,  N, 
and  the  average  single  scatter  properties  of  the  system  as 

IT 

g  *  1  -  <cos6>  ■  1  -  2x  /  p(0)  cosQ  sin0d6;  (1) 

o 

b  -  (NOs)-1  -  2*N  /irp(9)  sin©  d0;  (2) 

o 

t  -  Z/b,  (3) 

Z  being  the  geometrical  thickness  of  the  layer. 

The  scaling  relations  for  the  mean  and  standard  deviation  of 
the  path  length  according  to  Bucher^  are  then 

<i>  -  0.62  b  g0,94  T1,94  (4) 

and  ot  -  <( l  -  <t»2>  -  0.64  b  g0,81  t1,81  (5) 


while  the  beam  radius  at  which  the  intensity  drops  to  0.5  of  its 
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0.78  b  g 


-0.07 


t0.93 


and  the  fraction  of  the  incident  light  transmitted  through  the 
layer 


Ix/Io  *  1,69  (gT  +  1#42) 


(7) 


The  delay  in  the  arrival  time  of  an  incident  delta  function  pulse 
and  its  temporal  width  are  given  simply  by 
<At>  *  <*>/c 


_  n  -1  .  0.94  t1.94 

=  0.62  c  b  g  T 


and  o  =  o^/c 


**  -1  .  0.81  1.81 
0.64  c  b  g  x 


(8) 


(9) 


where  c  is  the  speed  of  light  in  the  medium. 

The  more  recent  analytic  and  computer  modelling  studies**”^ 
have  served  to  confirm  these  thick  cloud  results  and  provided  in 
addition  analytic  expressions  for  the  pulse  shape.  In  particular 
Ito  and  Furutsu^,  and  Ciervo^  have  expressed  the  temporal 
behavior  as  an  infinite  sum  of  decaying  exponentials  of  which  only 
the  first  few  terms  may  be  significant  under  most  circumstances. 
The  pulse  shape  has  also  been  modelled  as  either  a  gamma 
distribution  or  the  sum  of  two  gamma  distributions,^ 


I(t)  -  At  e"at  (10a) 

or  I(t)  *  Ait  e  ait  +  A2t  e  a2t .  (10b) 


Since  a  completely  detailed  comparison  of  pulse  shapes  is  beyond 
the  scope  of  the  present  work  only  the  gamma  distributions  will  be 
referred  to  below. 


III.  Scattering  System 
A.  Preparation 

The  experimental  measurements  reported  here  were  performed  on 
either  of  two  emulsions:  paraffin  oil  or  diiodomethane  (DIM)  in 
water.  The  first  was  chosen  for  its  ready  availability  and  ease 
of  preparation  and  handling;  the  second  for  its  index  of 
refract  ion ,  1.738,  which  has  nearly  the  same  refractive  index 
relative  to  water  as  water  does  to  air.  The  DIM/water  emulsion 
can  thus  closely  simulate  the  scattering  properties  of  a  cloud. 

In  each  case  lauryl  sulphate  was  used  as  a  surfactant  to  help 
stabilize  the  emulsion. 

All  materials  used  may  be  purchased  from  chemical  supply 
houses.  The  DIM  (CH2I2)  was  obtained  from  two  sources,  Alpha 
Products  and  Aldrich  Chemical  Co.  That  supplied  by  the  former  was 
very  clear  with  a  faint  yellowish  tinge  characterist ic  of  pure 
diiodomethane,  while  that  obtained  from  the  latter  had  been 
shipped  with  powdered  copper  added  for  stability  and  had  a  dark 
reddish  yellow  appearance  due  to  the  liberation  of  free  iodine. 

The  purity  of  the  DIM  was  restored  by  treating  it  with  an  aqueous 
solution  of  sodium  bisulphite  followed  by  anhydrous  sodium 
sulphate  to  remove  any  residual  water.  The  absorption  coefficient 
of  the  purified  DIM  was  measured  to  be  «  *  2.4  cm  1  at  a 
wavelength  of  532  ran.  The  lauryl  sulphate,  Na(CH2) 12HSO4 ,  was 
>  99%  GC  grade  manufactured  by  Fluka,  A.G.,  while  the 


paraffin  oil  was  purchased  from  VWR  Scientific,  Inc.  The  index  of 

refraction  of  the  oil  was  measured  to  be  1.4628  and  its  absorption 

coefficient  to  be  <  0.0004  cm  * .  All  water  used  including 

that  for  cleansing  was  deionized  by  a  reverse  osmosis  process  and 

18 

showed  a  bulk  resistivity  of  10  ohm-cm.  The  physical  properties 
of  the  paraffin  oil  and  diiodomethane  relevant  to  this  work  are 
listed  in  Table  I. 

The  emulsions  were  prepared  as  follows:  10  ml  of  oil  or  DIM 
in  1  i  of  0.5%  by  weight  solution  of  lauryl  sulphate  was  agitated 
vigorously,  circulated  through  a  stack  of  three  Nuclepore  filters 
(8  Un  pore  diameter)  for  several  minutes,  agglomerations  of  large 
drops  were  removed  from  the  surface  by  skimming  in  the  case  of  oil 
or  from  the  bottom  by  decanting  in  the  case  of  DIM,  This  process 
produced  a  coarse  emulsion  with  droplet  diameters  ranging  from  <  1 
Vm  to  50  Pm.  The  location  of  the  peak  of  the  size  distribution 
could  be  controlled  somewhat  by  varying  the  rate  of  flow  through 
the  filter  stack,  a  higher  flow  rate  shifting  the  peak  to  smaller 
diameters . 

The  coarse  emulsion  prepared  in  this  manner  could  be  stored 
for  periods  of  days  in  a  slowly  rotating  ("3  rpm)  cylindrical 
bottle  with  horizontally  oriented  axis  until  enough  had  been 
accumulated  for  the  next  stage  of  preparation. 

The  coarse  emulsions  were  further  refined  to  produce  more 
nearly  monodisperse  size  distributions  by  exploiting  the 
differential  drift  velocity  of  different  sized  drops  under  the 
influence  of  gravity.  According  to  Stokes  law  the  steady  state 


drift  velocity  of  spherical  body  in  a  viscous  medium  is  given  by 
V(d)  -  G  d2  (pD  -  p)/(18n)  (11) 

where  n  is  the  viscosity  of  the  medium,  G  is  the  acceleration  due 
to  gravity,  d  the  diameter  of  the  sphere,  and  p0  and  p  the 
density  of  the  medium  and  the  sphere  respectively.  With  water  as 
the  bulk  medium  h  ■  0.01  cp  and  p0  *  1.0  while  for  paraffin  oil 
p  *  0.8755  and  for  DIM  p  *  3.325.  Thus  a  1  Pm  diameter  oil 
droplet  should  rise  at  0.024  cm  hr  *,  and  a  10  Pm  diameter 
droplet  at  2.44  cm  hr  1 .  On  the  other  hand  a  1  Pm  diameter 
DIM  droplet  falls  at  0.46  cm  hr  1  and  a  10  Pm  DIM  droplet  at 
45.6  cm  hr  1 . 

A  separation  apparatus  consisting  of  a  20  cm  diameter 
plexiglas  cylinder  72  cm  long  and  a  hydraulically  activated 
"shutter"  located  8  cm  from  the  end  as  illustrated  in  Figure  1  was 
used  to  selectively  remove  either  large  or  small  droplets.  The 
cylinder  with  the  shutter  open  was  filled  with  coarse  emulsion 
prepared  according  to  the  procedure  described  earlier  and  placed 
in  a  temperature  controlled  (±  0.2°C)  mechanically  isolated 
cabinet  to  avoid  the  influence  of  convection  currents  and 
vibration.  After  a  time  calculated  to  allow  all  oil  droplets  of  a 
chosen  diameter,  dc,  to  rise  from  the  bottom  past  the  shutter, 

Tc  -  (72-8)/V(dc)  -  1152  n/[G|p0  -  p|dc2] ,  (12) 

the  shutter  was  closed.  The  procedure  followed  for  the  DIM 
systems  was  essentially  the  same  except  then  the  shutter  is 


located  near  the  bottom  of  the  column  and  the  droplets  fall  rather 
than  rise. 

The  short  portion  of  the  cylinder  contains  all  the  droplets 
of  the  size  chosen  and  greater  plus  the  droplets  of  smaller  size 
which  were  already  in  that  volume  or  which  were  able  to  move  past 
the  shutter  level  from  regions  near  it.  The  longer  part  of  the 
cylinder  contains  only  droplets  smaller  than  the  chosen  size.  A 
simple  calculation  using  Eqs.(ll)  and  (12)  and  the  cylinder 
dimensions  gives 

N.(d)  -  N  (d) [1  -  (d/d  )2]  ;  d  <  d 

L  O  C  C 

*  0  ;  d  >  d  (13) 

—  c 

N  (d)  -  1/9  N  (d) [1  +  8(d/d  )2]  ;  d  <  d 

s  o  c  c 

*  N  (d)  ;  d  >  d  (14) 

o  c  —  c 

where  NQ(d)  is  the  number  of  droplets  of  diameter  d  in  the 
coarse  emulsion  and  N^(d)  and  Ns(d)  are  respectively  the 
numbers  in  the  longer  and  shorter  lengths  of  the  column  after  time 
Tc. 

The  emulsion  in  either  length  of  the  cylinder  may  be 
collected,  made  up  to  its  original  volume  by  adding  water,  stirred 
to  redistribute  the  droplets  uniformly  throughout  the  volume  and 
the  process  repeated.  In  this  manner  it  was  possible  to 
significantly  reduce  the  variance  in  droplet  diameters.  For 
example,  assume  the  initial  size  distribution  is  uniform,  N0(d) 

*  N0,  and  oil  droplets  with  mean  diameter  near  10  Urn  are 


desired*  In  this  case  the  time  required  for  a  10  Pm  diameter 
droplet  to  rise  from  the  bottom  of  the  column  to  the  shutter  Tio  ■ 
26.23  hours.  If  the  emulsion  which  accumulates  in  the  long 
section  of  the  cylinder  is  discarded,  that  which  accumulates  in 
the  short  section  retained,  diluted  to  the  full  volume  of  the 
cylinder  and  the  process  repeated  say  four  times  the  resulting 
distribution  will  be 

Ni(d)  -  No9_4[ 1  +  8(d/10)2]4  ;  d  <  10 

-  N  ;  d  >  10  (15) 

o  — 

(d  in  micrometers).  If  then  one  final  stage  to  remove  the  drop¬ 
lets  with  diameter  greater  than  d  *  12  vm  is  performed  by  waiting 
for  a  time  T12  ■  18.21  hours  and  the  emulsion  in  the  long  section 
of  the  cylinder  retained,  the  size  distribution  which  results  is 

N2(d)  -  Ni(d)  [1  -  (d/12)2] 

-  N  9~4  [1  +  8(d/ 10) 2 ] 4  (1  -  (d/12)2]  ;  d  <  5 

o 

-  N  [1  -  (d/12)2]  ;  10  <  d  <  12 

o  — 

-  0  ;  d  >6  (16 

This  has  a  mean  diameter  <d>  ■  9.32  Pm  and  standard  deviation 
-  1.30  Pm. 

In  practice  several  factors  reduce  the  effectiveness  of  the 
separation  process.  There  is  a  tendency  of  particles  to  coalesce 
and  form  larger  droplets  although  this  is  much  reduced  by  the  use 
of  the  surfactant,  lauryl  sulphate.  Diiodomethane  is  slightly 


soluble  in  water  and  in  that  case  there  is  a  countervailing 
reduction  in  droplet  size  due  to  gradual  dissolution  of  the  DIM 
droplets.  This  effect  can  be  minimized  by  saturating  the  bulk 
medium  with  DIM.  Also  Stokes  law  applies  strictly  only  to  an 
isolated  sphere.  In  the  high  number  density  systems  of  interest 
the  movement  of  one  droplet  affects  those  around  it  so  that  some 
cooperative  motion  results.  On  the  other  hand  the  initial,  coarse 
emulsion  does  not  have  a  uniform  size  distribution  and  the  modal 
diameter  can  be  varied  to  advantage  by  means^f  the  flow  rate 
through  the  Nuclepore  filters.  Nevertheless,  it  is  possible  and 
practical  by  the  method  outlined,  to  produce  systems  of  droplets 
which  have  a  desired  mean  diameter  and  a  standard  deviation  on  the 
order  of  20%  of  the  mean.  In  the  case  of  oil  emulsions  a  period 
of  about  a  week  may  be  required  but  DIM  emulsions,  because  of  the 
greater  difference  in  density  and  consequent  greater  drift 
velocity,  may  be  prepared  in  the  course  of  a  day.  Oil  emulsion 
Code  F  (see  Table  II  and  Figure  2)  was  prepared  according  to  the 
detailed  program  outlined  above. 

B.  Characterization 

A  well  characterized  scattering  system  is  of  prime  importance 
if  the  experimental  results  are  to  be  of  more  than  limited 
usefulness.  The  preparation  technique  and  methods  of  handling  the 
emulsions  used  in  these  experiments  were  chosen  to  maintain  a 
stable  system  whose  characteristics  varied  sufficiently  slowly 


that  continuous  monitoring  of  the  system  was  unnecessary.  The 
emulsions  were  in  all  cases  thoroughly  mixed  and  kept  in  the 
rotating  scattering  cell  described  below  to  minimize  agglomeration 
of  the  droplets  and  concentration  gradients. 

Samples  of  the  emulsion  were  taken  at  intervals  throughout 
the  course  of  the  optical  measurements  and  the  size  distribution 
and  number  density  determined  by  direct  microscopic  examination. 
The  emulsion  sample  was  placed  in  a  cylindrical  observation  cell 
constructed  from  a  glass  microscope  slide,  a  0.25"  I.D.  by  0.01" 
stainless  steel  washer  and  a  glass  cover  slip.  After  the  droplets 
in  the  emulsion  had  risen  to  the  cover  slip  in  the  case  of  oil  or 
fallen  to  the  slide  in  the  case  of  DIM,  Polaroid  micrographs  were 
taken  at  400X.  At  this  magnification  and  resolution,  <  1  Pm,  the 
images  were  perfect  circles  indicating  that  the  droplets  were 
actually  spherical. 

The  diameter  of  each  individual  droplet  on  the  set  of 
micrographs  was  measured  directly,  classed  in  1  pm  intervals  and  a 
histogram  of  the  number  in  each  diameter  range  constructed.  The 
number  density  could  be  calculated  from  the  volume  represented  by 
each  micrograph  and  the  number  of  droplets  in  it.  In  all  cases, 
several  hundred  droplets  from  each  sample  were  measured  to  provide 
a  reasonable  statistical  average.  The  histogram  in  Figure  2 
was  constructed  in  this  manner  by  measuring  the  diameters  of  802 
droplets  of  oil  emulsion  Code  F.  This  had  mean  diameter  9.39  Pm, 
standard  deviation  1,69  Pm,  number  density  3  *  10  m  ,  and 
volume  fraction  0.14%.  Also  included  in  Figure  2  is  the 
corresponding  theoretical  size  distribution  given 


by  Eq.16.  Note  the  good  agreement  between  the  measured  and 
theoretical  mean  diameters  and  standard  deviations.  Table  II 
lists  these  quantitites,  the  number  density,  and  the  volume 
fraction  of  all  the  emulsions  used  in  the  optical  experiments. 

The  single  scattering  properties  of  the  emulsions  were 
determined  by  Mie  theory  from  the  size  distribution  and  bulk 
properties  of  the  paraffin  oil  and  diiodomethane •  The  scatter 
cross  section,  0g,  mean  of  the  cosine  of  the  scattering  angle, 
<co86>,  and  the  albedo,  were  calculated  for  droplets  with 
diameters  ranging  from  0.1  Pm  to  30  ttn  and  the  average  of  these 
quantities  weighted  by  the  size  distribution  histogram  determined 
for  each  emulsion.  The  scattering  properties  of  each  emulsion  are 
presented  in  Table  III. 

C.  Scattering  Cell 

The  scattering  cell  is  illustrated  in  Figure  3.  A  40  cm 
diameter  by  12  cm  cylindrical  aluminum  tank  is  suspended  in 
bearings  with  its  axis  horizontal  and  provision  made  for  it  to  be 
rotated  at  a  constant  speed  of  about  1  rpra.  One  end  of  the  tank 
is  closed  with  a  plate  glass  exit  window,  the  other  with  an 
aluminum  plate  with  a  5  cm  circular  hole  in  its  center  through 
which  a  sliding  hollow  tube  protrudes.  The  end  of  the  tube  inside 
the  tank  is  closed  with  a  2.5  cm  diameter  glass  entrance  window 
and  is  attached  to  a  38  cm  diameter  circular  disc  which  acts  as  a 
false  tank  end.  The  entire  surface  of  the  tank,  apart  from  the 
windows,  is  painted  flat  black  to  avoid  light  being  reflected  back 
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into  the  scattering  medium  once  it  has  'escaped.'  It  is  thus 
possible  by  varying  the  position  of  the  entrance  window  and 
associated  false  end  relative  to  the  exit  window  to  study 
scattering  in  a  medium  bounded  by  two  "infinite"  parallel 
absorbing  planes  separated  by  a  distance  ranging  from  2  to  117  mm. 
D.  Optical  System 

The  source  of  the  optical  pulses  for  the  delay  and  stretching 
measurements  was  a  passively  mode-locked  Nd:YAG  laser  which 
produced  a  train  of  around  20  pulses,  each  35  ps  in  duration  and 
carrying  *  100  uJ  energy.  The  1064  nm  radiation  was  double 
passed  through  an  angle  tuned  KDP  crystal  converting  approximately 
20Z  of  the  energy  to  second  harmonic  532  nm  light  which  was 
separated  from  the  infrared  fundamental  by  a  dichroic  mirror  (see 
Figure  4).  A  single  green  pulse  was  selected  from  the  train  by 
means  of  a  Pockels  cell,  the  selected  pulse  being  ~  25  ps  in 
duration. 

A  portion  of  the  energy  of  the  selected  pulse  was  split  off 
by  a  beam  splitter  to  trigger  the  detection  system  and  a  further 
portion  split  off  and  routed  through  a  delay  line,  around  the 
scattering  cell,  to  the  detector  to  provide  an  accurately  timed 
reference  event.  The  major  portion  of  the  selected  pulse  was 
directed  onto  the  entrance  window  of  the  scattering  tank. 

Provision  was  made  to  place  neutral  density  filters  in  the  path  of 
either  the  reference  or  primary  beam  to  allow  both  the  reference 
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and  the  scattered  pulses  to  be  recorded  within  the  dynamic  range 
of  the  detection  system*  For  all  measurements  except  those  on 
emulsion  Code  A  a  Cornu  pseudodepolarizer  was  placed  in  the  beam 
to  provide  unpolarized  incident  pulses* 

The  detector  employed  in  all  the  pulse  measurements  was  a 
Hamamatsu  Model  C979  streak  camera  capable  of  resolving  events 
separated  in  time  by  less  than  10  picoseconds  when  operated  at  its 
fastest  streak  speed  and  at  lower  streak  speeds  it  can  display 
event 8  somewhat  longer  than  5  nanoseconds.  The  angular  acceptance 
of  the  streak  camera  is  approximately  cosine  squared  with  a  36* 
full  angle  at  half-power.  One  series  of  measurements  on  emulsion 
Code  G  was  made  with  an  aperture  restricting  the  field  of  view  to 
16*. 

The  streak  camera  provides  both  analog  and  digital  readout 
with  the  actual  streak  and  its  intensity  profile  being  displayed 
on  a  television  monitor.  The  immediately  available  record  is 
invaluable  in  setting  up  and  aligning  the  experiments.  The 
digital  data  record  is  readily  stored  on  magnetic  tape  for 
subsequent  numerical  analysis. 

Some  spatial  beam  broadening  measurements  were  made  by 
directing  the  beam  from  a  cw  argon  ion  laser  operating  at  514  nm 
onto  the  entrance  window  of  the  scattering  tank.  Beam  profiles  at 
the  exit  window  were  obtained  by  scanning  a  pin  diode  detector 
along  a  diameter  of  the  exit  window. 


IV.  Experimental  Results 

A.  Pulse  Stretching,  Delay  and  Attenuation 

The  pulse  scattering  experiments  were  all  performed  according 
to  the  following  procedure.  After  the  emulsion  was  prepared  and 
characterised  the  scattering  cell  was  filled  and  placed  in  the 
beam  as  illustrated  in  Figure  4.  The  movable  entrance  window  to 
the  cell  was  positioned  at  the  desired  distance  from  the  exit 
window  and  a  set  of  approximately  20  intensity  versus  time  records 
of  the  reference  and  scattered  pulses  were  obtained  with  the 
streak  camera.  The  digitized  events  were  recorded  on  magnetic 
tape  for  subsequent  computer  analysis.  The  entrance  window  was 
then  repositioned  to  a  different  path  length  and  the  procedure 
repeated. 

The  digitized  records  of  the  individual  events  for  each  path 
length  were  put  in  precise  register  by  taking  the  midpoint  of  the 
half  maximum  line  of  each  reference  pulse  as  time  zero  and  summing 
the  20  individual  records.  This  minimized  the  effects  of  pulse  to 
pulse  variation  and  reduced  random  errors.  Figures  5  and  6  are 
typical  examples  of  these  composite  records  illustrating  the 
dramatic  pulse  stretching  phenomena  due  to  multiple  scattering  in 
optically  thick  media.  Figure  5,  the  results  obtained  for  an 
optical  thickness  t  ■  9,  shows  the  scattered  pulse  to  be  only 
slightly  broader  than  the  reference  pulse  while  Figure  6,  optical 
thickness  t  ■  63,  shows  considerable  stretching  of  the  scattered 


pulse  relative  to  the  reference  pulse. 

The  integrated  intensity  (area),  mean  arrival  time 
(centroid),  and  rms  pulse  width  of  both  the  scattered  and 
reference  pulses  were  calculated  from  the  composite  records  for 
each  experimental  condition.  From  the  integrated  intensities, 
knowledge  of  the  neutral  density  filters  in  each  path  and  the 
corresponding  intensities  measured  in  the  absence  of  scattering, 
i.e.,  with  pure  water  in  the  scattering  cell,  the  fraction  of  the 
incident  intensity  which  was  collected  by  the  receiver  could  be 
calculated. 

In  a  similar  manner  the  displacement  of  the  centroid  of  the 
scattered  pulse  from  that  of  the  reference  could  be  used  to 
determine  the  excess  delay  due  to  scattering.  This  entailed 
taking  into  account  the  optical  elements,  e.g.  neutral  density 
filters,  in  each  path,  the  extra  distance  the  scattered  pulse  had 
to  travel  through  water  at  the  longer  paths  and  the  amount  of 
delay  introduced  deliberately  in  the  delay  line  to  position  the 
two  pulses  a  convenient  distance  apart  on  the  streak  record. 

It  may  be  noted  that  the  delay  measurements  reported  here  are  the 
first  ever  measured.  Prior  experiment s^ *3,4,6  have  not  been 
able  to  determine  delays  precisely  because  of  the  lack  of  an 
adequately  timed  reference  signal. 

The  pulse  stretching  caused  by  multiple  scattering  was 
calculated  by  subtracting  the  second  moment  of  the  reference  pulse 
from  that  of  the  scattered  pulse,  taking  the  square  root  and 
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doubling.  Twice  the  rms  width,  2  crt,  was  chosen  for  comparison 
purposes  since  it  is  closer  to  the  7WHM  of  a  pulse.  For  example 
the  FWHM  of  a  gauss ian  pulse  is  2.35  o ^  while  that  of  a  te~a* 
shaped  pulse,  is  1.73  0t. 

The  delay  time,  At;  pulse  stretching,  2  ot;  and  the  ratio 

of  the  received  intensity  to  the  incident  intensity  determined  in 

the  manner  described  for  each  path  length  (optical  thickness)  and 

for  all  of  the  emulsions  listed  in  Tables  II  and  III  are  presented 

in  Tables  IVA  through  IVG.  Also,  the  delay  time  results  are 

plotted  in  Figures  7  and  8,  the  pulse  stretching  results  in 

Figures  9  and  10,  and  the  relative  intensity  results  in  Figures  11 

and  12.  In  order  to  facilitate  intercomparison  of  the  results  the 

delay  and  pulse  stretching  measurements  have  been  normalized  by 

multiplying  them  by  c  b  1  g  ^#^Vo.62  and  c  b  1 
*”0.81 

g  *  / 1 . 28  respectively  with  the  values  of  b  and  g  used 

being  those  appropriate  to  the  particular  emulsion. 

Theoretically,  according  to  Eqs.(8)  and  (9),  the  data  should  then 
fall  on  the  lines  t  *  in  the  case  of  the  delay  measurements 
and  on  the  line  x  '  for  the  pulse  stretching  measurements. 

Again  for  the  purpose  of  ease  of  comparison  the  relative 
intensities  are  plotted  against  gT  since  according  to  Eq.(7)  the 
points  should  lie  along  1.69/(gx  +  1.42). 

B.  Field  of  View  Effect 

A  series  of  pulse  scattering  measurements  was  also  made  with 
the  field  of  view  of  the  streak  camera  restricted  to  16*  full 
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angle  rather  than  its  unrestricted  field  of  view  of  36*.  These 
measurements  were  made  in  the  same  way  as  described  above  but  on 
only  one  emulsion,  Code  G.  The  delay,  pulse  stretching  and 
relative  intensity  determined  for  four  optical  thicknesses  are 
displayed  in  Table  V  and  are  also  plotted  as  the  triangular  points 
in  Figures  8,  10  and  12  respectively.  These  should  be  compared 
directly  with  those  results  for  the  same  optical  thicknesses  in 
Table  IVG  and  the  corresponding  diamond  shaped  points  plotted  in 
Figures  8,  10,  and  12. 

C.  Beam  Spreading 

The  irradiance  profile  at  the  exit  window  of  the  scattering 
cell  was  mesured  by  scanning  a  wide  field  of  view  (>  120*)  pin 
photodiode  across  the  front  of  the  cell.  The  profiles  obtained 
when  the  entrance  window  was  1,  4,  6,  8  and  11.7  cm  from  the  exit 
window  and  the  tank  was  filled  with  emulsion  Code  G  are  plotted  in 
Figure  13.  The  measured  beam  radii  (half  power  point)  for  these 
curves  are  listed  in  Table  VI  and  plotted  in  Figure  14  along  with 
the  theoretical  line  rc  *  0. 78  b  g  0 ’ 7  t° ’ 93  from 
Eq.(6). 
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V.  Discussion 

A  casual  examination  of  Figure  7  shows  quite  reasonable 
agreement  between  the  theoretical  values  of  the  delay  and  those 
measured  in  the  scattering  experiments  on  the  diiodome thane/ water 
emulsions  Codes  B,  C  and  D.  There  is  some  discrepancy  at  optical 
thicknesses  less  than  about  20,  but  at  small  optical  thickness 
the  probable  error  is  of  the  order  of  -  10  ps  and  as  illustrated 
by  the  error  bars  the  disagreement  is  not  excessive.  All  the  data 
from  emulsion  Code  A,  however,  falls  somewhat  below  and  parallel 
to  the  1.94  power  law  line,  a  fact  which  has  defied  rational 
explanation.  These  measurements  were  made  with  linearly  polarized 
incident  pulses  but  the  detection  system  was  not  polarization 
sensitive  and  no  effect  should  result. 

The  pulse  stretching  data  for  the  diiodomethane  emulsions 
displayed  in  Figure  9  shows  behavior  similar  to  the  corresponding 
delay  measurements.  The  increased  scatter  is  due  to  the 
inherently  poorer  accuracy  of  the  second  moment  calculations  and 
the  greater  weight  accorded  to  the  tail  of  the  scattered  pulse. 

The  fall  off  of  the  data  points  below  the  theoretical  line  at 
large  optical  thickness  may  be  due  to  two  factors.  First,  the 
duration  of  the  scattered  pulses  is  in  some  cases  so  great  that 
the  entire  event  could  not  be  captured  within  the  time-scale  range 
of  the  streak  camera.  Those  events  in  which  the  tail  of  the 
scattered  pulse  is  known  to  be  lost  are  identified  in  the  tables 
by  a  dagger.  Second,  the  single  scatter  albedo  calculated  for  the 
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diiodomethane  emulsions  is  less  than  0.9997  and  at  long  path, 

large  optical  thickness  absorption  of  the  scattered  light  may 

reduce  the  measured  pulse  stretching  and  delay. 

On  the  other  hand  there  is  poor  agreement  between  the 

theoretical  delay  and  stretching  and  the  measured  values  for  the 

paraffin  oil/water  emulsions  (Codes  E,  F,  and  G)  as  illustrated  in 

Figures  8  and  10.  The  experimental  data  all  falls  considerably 

below  the  theoretical  lines  and  moreover  the  functional  dependence 

on  optical  thickness  does  not  appear  to  be  a  simple  power  law.  m 

particular  the  experimental  delay  data  seem  to  indicate  a  very 

slow  initial  increase  with  optical  thickness  on  the  order  of 
0*3 

t  *  for  t  <  20  and  a  faster  rise  approximately  proportional  to 
2*6 

T  for  t  >  30.  The  one  scattering  parameter  which  is 
significantly  different  for  the  two  types  of  emulsion  is  the 
asymmetry  factor  g  =  1  -  <cos0>  which  is  as  large  as  0.17140  for 
the  DIM  emulsions  and  as  small  as  0.05045  for  the  oil  emulsions. 

It  is  likely  that  the  simple  dependence  on  g  assumed  here  is 
inadequate  and  that  higher  moments  of  cos  6  must  be  included  to 
describe  multiple  scattering  in  systems  with  very  strong  forward 
scattering,  g  +  0.  It  is  worth  noting  that  there  is  more 
satisfactory  agreement  for  the  diiodomethane  emulsions  which 
exhibit  the  same  degree  of  foward  scattering  as  water  droplets. 

The  effect  on  the  delay  and  stretching  caused  by  a  reduced 
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field  of  view  is  illustrated  in  Figures  8  and  10  by  the  triangular 
data  points.  These  measurements  were  made  on  emulsion  Code  G  and 
should  therefore  be  compared  most  closely  with  the  diamond  shaped 
data  points  at  the  same  optical  thickness.  A  dramatic  decrease  in 
the  measured  delay  and  stretching  is  evident.  The  delay  is 
reduced  by  a  factor  of  4  and  the  stretching  by  a  factor  of  5  at  an 
optical  thickness  of  73. 

The  variation  of  the  integrated  intensity  of  the  scattered 
radiation  with  optical  thickness  displayed  in  Figures  11  and  12 
does  show  the  functional  form  predicted, 

IT/I0  *  (gT  +  1.42)  1 .  However,  the  constant  of 

proportionality  is  10  times  smaller  than  expected  if  the  detector 
were  located  at  the  surface  of  the  exit  window.  In  these 
experiments  the  streak  camera  entrance  slit  was  located  20  cm  from 
the  glass  exit  window  a  fact  which  may  account  in  part  for  this 
difference. ^  The  integrated  intensity  measurements  were  the 
least  accurate  of  any  made  during  these  experiments  but  are 
certainly  correct  to  within  a  factor  of  2,  i.e.,  ±  0.3  on  the 
logarithmic  scale. 

The  effect  of  the  reduced  field  of  view  on  the  received 
intensity  is  clear  (triangles  in  Figure  12).  The  data  falls  off 
more  rapidly  with  increasing  optical  thickness  as  would  be 
expected  from  the  fact  that  less  of  the  light  scattered  through 
extreme  angles  is  accepted.  The  down  turn  in  the  intensity  data 


for  emulsion  Code  A  evident  at  very  large  optical  thicknesses  in 
Figure  11  is  probably  another  manifestation  of  the  effects  of  the 
reduced  albedo  of  this  large  droplet  diiodomethane  emulsion. 

The  irradiance  profiles  shown  in  Figure  13  are  similar  to  the 
predicted  shape. 2  The  measured  half  power  beam  radii ,  Table  VI 
and  Figure  14,  are  however  much  smaller  than  the  0.78  b  g 
t  *  predicted,  especially  at  small  optical  thickness.  These 
measurements  were  made  on  the  paraffin  oil/water  emulsion  Code  G 
and  the  discrepancy  can  again  be  attributed  to  the  extremely 
strong  forward  scattering,  g  ■  0.05045,  which  results  from 
droplets  with  large  size  parameters  and  small  relative  index  of 
refraction. 

A  perfunctory  study  of  the  shape  of  the  scattered  pulses  has 
been  made.  First  an  attempt  was  made  to  fit  the  functional  form 
of  a  gamma  distribution,  Eq.(lOa),  to  the  scattered  pulse  shape. 
The  goodness  of  fit  was  so  poor  however  that  this  was  abandoned 
and  the  two  gamma  distribution  form,  Eq.(lOb),  currently  in 
vogue, ®  was  tried.  This  too  yielded  poor  fits  to  pulses 
scatterd  by  optically  thin  systems  and  only  moderately  good  fits 
to  long  path,  optically  thick  data  as  illustrated  in  Figures  15 
and  16.  These  are  least  squares  fits  of  the  2  gamma  form  to  data 
obtained  from  the  experiments  on  the  diiodomethane  emulsion  Code  A 


and  the  paraffin  oil  emulsion  Code  G  at  optical  thicknesses  T  ■ 

167  and  x  *  107  respectively. 

The  relatively  poor  agreement  between  the  theoretical  and 
observed  pulse  shapes  indicates  that  not  much  is  to  be  gained  by 
using  the  fitted  curve  to  calculate  the  delay  or  stretching  except 
in  those  cases  where  the  tail  of  the  scattered  pulse  extended 
beyond  the  range  of  the  streak  camera  time  scale.  The  equation  of 
the  theoretical  curve  in  Figure  15  is  I(t)  ■  0.94  t  exp  (-1.94  t) 

+  0.27  t  exp(-0.60  t),  t  in  nanoseconds.  Twice  the  rms  width  of 
of  this  pulse,  2  *  4.60  ns.  This  is  to  be  compared  to  the 

directly  calculated  value  of  1.74  ns  (see  Table  IV  A).  When  this 
is  normalized  by  multiplying  by  c  b  1  g  0#81/1*28,  with 
b  and  g  appropriate  to  emulsion  Code  A,  one  gets  7800  for  the 
normalized  stretching  which  when  plotted  as  the  solid  square  datum 
on  Figure  9  is  still  somewhat  below  the  theoretical  line  for  that 
optical  thickness.  It  is  however  much  closer  than  the 
corresponding  hollow  square. 

The  principal  conclusion  which  can  be  drawn  from  this  work  is 
that  there  is  reasonable  agreement  between  the  Monte  Carlo 
simulation  results  and  these  experiments  for  the 
diiodomethane/water  systems.  Since  these  systems  mimic  water 
droplet  clouds  the  application  of  the  theory  to  scattering  in 
clouds  and  fogs  is  certainly  in  order.  The  paraffin  oil/water 
results  on  the  other  hand  indicate  that  the  present  theory  is 
inadequate  for  systems  exhibiting  extreme  forward  scatter  and  a 


more  complete  theory  is  needed  perhaps  involving  higher  moments? 

o 

of  the  scattering  function,  <cos  6>,  etc. 

There  is  of  course  much  work  which  must  be  done  before  the 
whole  story  of  the  effects  of  multiple  scattering  on  optical  pulse 
propagation  is  known.  A  wider  range  of  index  of  refraction  should 
be  explored  as  well  as  an  expanded  range  of  size  parameters.  The 
effects  of  the  field  of  view  of  the  receiver  needs  to  be 
thoroughly  studied  as  does  the  temporal  behavior  of  polarized  and 
depolarized  scattered  radiation  when  the  incident  pulse  is 
linearly  polarized.  A  complete  pulse  shape  study  would  also 
provide  a  valuable  addition  to  the  understanding  of  multiple 
scattering. 
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TABLE  I 


Paraffin  Oil 

Diiodomethane 

Density 

0.8755 

3.325 

Index  of  Refraction 

1.4628 

1.738 

Index  Relative  to  Water 

1.0974 

1.304 

Absorption  Coefficient 

<  4  x  io”4  cm”1 

2.4  cm"1 

@  532  nm 

Drift  Velocity  of 

2.44  cm  hr  1 

-45.6  cm  hr  1 

10  Vtn  Droplet 

in  water 

TABLE  Ii 

EMULSION  PROPERTIES 

Number  Density 
(m"3) 


Emulsion  Mean  Diameter  Standard 
Code  and  Type  (Pm)  Deviation  (Pm) 


Volume  Fraction 
(percent) 


TABLE  III 


SCATTERING  PARAMETERS 


Emulsion 

Code 

Size 

Parameter 
*<d>/ X 

Mean  Scatter 
Cross  Section 
o#  (um)Z 

Scattering 
Length 
b  (mm) 

Asymmetry 

Factor 

g  ■  1  -  <coa8> 

Single 

Scatter 

Albedo 

A 

* 

65.1 

148 

1.34 

0.14315 

0.99926 

B 

26.6 

24.2 

0.389 

0.17073 

0.99970 

C 

30.5 

36.0 

3.40 

0.17140 

0.99964 

D 

62.5 

130 

0.453 

0.14177 

0.99926 

E 

14.6 

8.95 

1.27 

0.07788 

1.00000 

F 

73.9 

152 

2.19 

0.06075 

1.00000 

G 

101.8 

298 

1.09 

0.05045 

1.00000 

•ft 

Wavelength  in  medium  (water)  X  *  0.532/1.333  *  0.399  ym 


TABLE  IVA 


Path  Length  Optical  Delay  Stretching  Relative 

Z  (cm)  Thickness  At  (ps)  2  (ps)  Intensity 


TABLE  IVB 


Path  Length 
Z  (cm) 


Optical 

Thickness 


Delay 
At  (  ps  ) 


Stretching 

2 


Relative 

Intensity 


TABLE  IVE 


Path  Length 
Z  (cm) 


Opt ical 
Thickness 


Delay 
At  ( ps ) 


Stretching 

2  a 


Relat ive 
Intensity 


(‘ 


TABLE  IVF 


Path  Length 

Z  (cm) 

Optical 

Thickness 

T 

Delay 

At  (ps) 

Stretching 

2  °t 

(p8) 

Relative 

Intensity 

0.5 

2 

13 

17 

1.5  x  10_1 

1.0 

5 

2 

18 

3.6  x  IQ-2 

2.0 

9 

11 

17 

CO 

1 

O 

—4 

X 

00 

«— 4 

4.0 

18 

34 

17 

8.9  x  10"5 

6.0 

27 

33 

35 

4.5  x  io'5 

8.0 

37 

64 

79 

3.1  x  io-5 

10.0 

46 

95 

130 

1.9  x  i(f5 

11.7 

53 

120 

170 

1.1  x  io-5 

38 


k  * 

l 

l 
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Figure  1.  Apparatus  for  separating  large  and  small  buoyant 
droplets  in,  e.g.,  an  oil/water  emulsion.  A  thin 
stainless  steel  shutter  may  be  drawn  across  the  column 
to  isolate  the  top  1/9  of  the  volume  from  the  bottom 
8/9.  For  diiodoraethane/water  or  other  emulsions  with 

r4 

r  droplet 8  more  dense  than  the  bulk  medium  the  column  is 


inverted  with  the  shutter  located  near  the  bottom. 


DROPLET  DIAMETER  l/im) 

Figure  2.  Size  distribution  histogram  for  emulsion  Code  F.  The 

dashed  line  is  the  theoretical  distribution  according  to 
Eq.  16.  The  measured  mean  diameter  and  standard 
deviation  are  9.39  Un  and  1.69  vm  compared  to  the 
predicted  values  of  9.32  Um  and  1.30  wm  respectively.  A 
total  of  802  droplet  diameters  were  measured  on  400X 
micrographs  of  samples  of  the  emulsion. 
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Window 


Figure  3-  Scattering  cell.  The  cylinder  rotates  at  a  constant 
speed  of  1  rpm  to  prevent  the  emulsion  droplets  from 
collecting  at  the  top  or  bottom.  The  distance  between 
the  entrance  and  exit  windows  is  variable  between  0.2  cm 
and  11.7  cm.  All  interior  surfaces  are  flat  black  to 
reduce  reflection  of  escaping  radiation  back  into  the 
scattering  medium. 
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Figure  5,  Typical  intensity  versus  time  streak  camera  record.  The 
left  hand  pulse  is  the  reference  and  the  right  hand  the 
pulse  which  has  passed  through  the  scattering  medium. 
Optical  thickness  t  *  9,  emulsion  Code  A. 


f/ 
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Figure  6.  The  same  as  Figure  5  except  t  *  63. 
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Figure  7.  Delay  in  the  mean  arrival  time  due  to  multiple 

scattering  versus  optical  thickness  for  diiodomethane 
emulsions.  Measured  delay  times  have  been  normalized 
by  multiplying  by  c  b  1  g  °'^4/0.62.  □  -  code  A;  0  - 
code  B;  0-  code  C;  and  A  -  code  D.  The  line  is 
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NORMALIZED  PULSE  WIDTH 


restricted  field  of  view  (16°)  detector. 


Figure  11.  Logarithm  of  the  relative  integrated  intensity  of 

scattered  pulses  versus  optical  thickness  times  the 
asymmetry  factor,  gT,  for  diiodomethane  emulsions. 

□  -  code  A;  0  -  code  B;  0  -  code  C;  and  A  -  code  D. 
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O  12345678 

gT 

Figure  12.  The  same  as  Figure  11  but  for  paraffin  oil  emulsions. 

o  -  code  E;  0  -  code  F;  0  -  code  G;  and  A  -  code  G 
with  restricted  field  of  view  (16*)  detector. 
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OPTICAL  THICKNESS,  t 

Beam  radius  versus  optical  thickness  as  measured  from 
the  half  power  points  on  the  profiles  displayed  in 
Figure  13.  The  line  is  r  *  0.78  b  g  0-7  t0,93. 


Figure  15.  Least  squares  fit  of  two  gamma  distribution  functional 
form  to  pulse  scattered  by  167  optical  thicknesses  of 
diiodomethane  emulsion  code  A.  -  0.94  t  exp(-1.94 


t)  +  0.27  t  exp(-0.60  t),  t  in  nanoseconds. 


Figure  16.  The  same  as  Figure  15  but  for  a  pulse  scattered  by  107 
optical  thicknesses  of  paraffin  oil  emulsion  code  G. 
-  3.28  t  exp(-2.78  t)  ♦  0.34  t  exp(-0.77  t),  t  in 


nanoseconds. 


